Averaging of the electron effective mass in multicomponent transparent conducting 

oxides 
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We find that layered materials composed of various oxides of cations with s 2 electronic config- 
uration, XI2O4, X=In or Sc, Y—G&, Zn, Al, Cd and/or Mg, exhibit isotropic electron effective 
mass which can be obtained via averaging over those of the corresponding single-cation oxide con- 
stituents. This effect is due to a hybrid nature of the conduction band formed from the s-states of all 
cations and the oxygen p-states. Moreover, the observed insensitivity of the electron effective mass 
to the oxygen coordination and to the distortions in the cation-oxygen chains suggests that similar 
behavior can be expected in technologically important amorphous state. These findings significantly 
broaden the range of materials as efficient transparent conductor hosts. 



Transparent conducting oxides (TCO) - the vital part 
of optoelectronic devices - have been known for a century 
and employed technologically for decades [l|, 0, H, 0, 
Yet, the current TCO market is dominated by only three 
materials, ln 2 03, Sn02 and ZnO, and the research efforts 
are primarily focused on the oxides of post-transition 



metals with (n — l)d w ns 2 electronic configuration. De- 
spite excellent optical and thermal properties as well as 
low cost, oxides of the main group metals, such as AI2O3, 
Si02, MgO and CaO, have never been considered as can- 
didates to achieve useful electrical conductivity due to 
the challenges of efficient carrier generation in these wide- 
bandgap materials [1, 0, HI • 

Multicomponent TCO with layered structure, e.g., 
InGaZn0 4 % EH, ED, El ES El , drew attention due to a 




GaZnO 



possibility to separate carrier donors (traditionally, oxy- 
gen vacancies or aliovalent substitutional dopants) and 
the conducting layers where carriers are transfered with- 
out charge scattering on the impurities. In lnGaZn04, 
octahedrally coordinated In layers alternate with double 
layers of oxygen tetrahedrons around Ga and Zn, Fig. Q] 
Because octahedral oxygen coordination of cations was 
long believed to be essential for a good transparent con- 
ductor d, El 03, 

53, EH , it has been suggested that in 
lnGaZn04 the charge is transfered within the InOi.5 lay- 
ers while the atoms in GaZn02 5 layers were proposed as 



; 5 layers w 

candidates for efficient doping [§, E3, EH 
has been argued that InGaOa(ZnO) 
tor El- 



Conversely, it 
is a Zn 4s conduc- 



To understand the role of local symmetry in the intrin- 
sic transport properties of TCO's and to determine the 
functionality of structurally and chemically distinct lay- 
ers in InGaZn04, we employ ab-initio density functional 
approach to study the electronic properties of various sin- 
gle and multi-cation oxides. Further, using InGaZn04 as 
a test model which examplifies not only the structural 
but also combinatorial peculiarities of complex TCO's, 
we survey other ns 2 cations - beyond the traditional In, 
Zn and Ga - for a possibility of being effectively incor- 
porated into novel multicomponent TCO hosts. 

Isotropy of the electronic properties in InGaZnO^. 



FIG. 1: The unit cell of InGaZnCU where three similar blocks 
consisting of one InOi.5 and two GaZn02.5 alternate along the 
[0001] direction. Ga and Zn atoms are distributed randomly. 
Three-dimensional interatomic (background) charge density 
distribution is evident from the contour plot calculated in the 
(Oil) plane. The plotted charge density corresponds to the 
carrier concentration of ~lxl0 18 cm~ 3 . 



The electronic band structure calculations for InGaZn04 
show that the atoms from both InOi.5 and GaZn0 2 .5 
layers give comparable contributions to the conduction 
band, fig. [2j leading to a three-dimensional distribution 
of the charge density, fig. [T] Moreover, the isotropy of 
the electronic properties in this layered material mani- 
fests itself in the electron effective masses being nearly 
the same in all crystallographic directions (table HI . 

The conduction band in InGaZn04 consists of a set 
of highly dispersed parabolic bands, fig. [21a). Since the 
band gap values in the corresponding single metal oxides 
are different, one may expect that each band is attributed 
to a certain cation in this multicomponent compound. 
However, we find that each band cannot be assigned to 
a state of a particular atom since all atoms in the cell, 
including the oxygen atoms, give non-negligible contri- 
butions to the conduction band wavefunction, fig. [2] and 
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FIG. 2: Partial density of states at the bottom of the conduc- 
tion band for InGaZnCU and InAIMgO^. Atoms from both 
In-Ol and Ga-Zn-02 (or Al-Mg-02) layers give non-negligible 
contributions. 



table |H 

The conduction band dispersion calculated along the 
[0001] crystallographic direction for a single unit cell, 
fig- Etb) , reveals that the multiple bands can be at- 
tributed to a "folding" of one parent band. Triple unfold- 
ing of the conduction band corresponds to the three-time 
reduction (expansion) of the conventional unit cell (Bril- 
louin zone) in the z direction. Since the block of three 
layers (one InOi.5 and two GaZnC>2.5 layers, fig. [TJ is 
repeated in the unit cell via translation, the splitting be- 
tween the resulting bands [for the fc-vector equal to — , — , 
— , fig. EJc)] is negligible. Although the subsequent 
unfolding into three individual layers is not justified be- 
cause the three layers are structurally and chemically dis- 
similar, fig. [1] we find that the band can be unfolded 
again, fig. [3{c) . The resulting highly dispersed band is in 
accord with the stepless increase of the density of states, 
fig- Eld). Thus, the conduction band can be "unfolded" 
nine times that corresponds to the total number of layers 
in the unit cell. Therefore, the electronic properties of 
the individual layers are similar. 

Unconventional s 2 -cations at work. A two-dimensional 
electronic structure could be expected in InAlMg04 and 
ScGaZnC^ since the band gap values in SC2O3, AI2O3 
and MgO are at least twice larger than those in In203, 
Ga2C>3, CdO and ZnO and, hence, the unoccupied s- 
states of Sc, Al and Mg should be located deeper in the 
conduction band. From the analysis of the partial den- 
sity of states for InGaMgC>4, InAlCd04, ScGaZn04 and 
InAlMg04, we find that although the contributions to 
the bottom of the conduction band from Sc, Al and Mg 
atoms are notably reduced, these states are available for 
electron transport, fig. El Thus, similar to InGaZn04 
where the cations s-states are energetically compatible, 
in all multicomponent oxides considered, the conduction 



band wavefunction is a combination of the s-states of all 
cations and the p-states of the oxygen atoms. The con- 
tributions from the chemically distinct layers are compa- 
rable, table |T1 and, consequently, these complex oxides 
exhibit three-dimensional network for the electron trans- 
port and isotropic electron effective mass, table HI 

Comparison to single-cation TCO's. The unfolded 
conduction band in the layered multicomponent mate- 
rials resembles those of single-cation TCO's, e.g., In 2 03, 
cf. Figs. [3] (c) and (e). Such a highly-dispersed single 
conduction band is the key attribute of any conventional 
n-type TCO host d M, 0, ED, lH ■ Upon proper 
doping, it provides both high mobility of extra carriers 
(electrons) due to their small effective mass, and low 
optical absorption due to high-energy inter-band tran- 
sitions from the valence band, E^, and from the par- 
tially filled conduction band, E c , fig. [4] Even in rela- 
tively small bandgap oxides, e.g., CdO where the optical 
band gap, E g , is 2.3 eV, the high energy dispersion en- 
sures a pronounced Fermi energy displacement with dop- 
ing (so called Burstein-Moss shift) which helps to keep 
the intense transitions from the valence band out of the 
visible range. However, large carrier concentrations re- 
quired for good electrical conductivity may result in an 
increase of the optical absorption due to low-energy tran- 
sitions from the Fermi level up into the conduction band 
as well as plasma frequency. Application-specific optical 
properties and desired band offsets (work functions) can 
be attained in a multicomponent transparent conductor 
with a proper composition. 

In both single and multi-cation oxides, the conduc- 
tion band is formed from the empty s-states of the metal 
atoms and the oxygen antibonding p-states, e.g., fig. [2] 
For multicomponent oxides we find that even at the bot- 
tom of the conduction band, i.e., at T point, the contri- 
butions from the oxygen p-states are significant, table |U 
Thus, the key feature of the conduction band in a conven- 
tional TCO - its high energy dispersion - originates in 
a strong interaction between the cation s-states and the 
anion antibonding p-states 17J. The direct s-s overlap 



is insignificant, fig. [TJ and therefore, the s-s interactions 
which were commonly assumed to play to key role in the 
electronic properties of TCO [H, 0, do not gov- 
ern the transport properties in these oxides. Indeed, the 
small electron effective mass in s 2 -cation oxides is deter- 
mined by the strong s-p interactions [26j [. 

From orbital symmetry considerations, one can see 
that the oxygen coordination of cations does not affect 
the s-p overlap. Instead, the largest overlap should be at- 
tained in materials where the oxygen atom is coordinated 
octahedrally by cations with the extended s-orbitals [17 1. 
Our systematic comparison of the calculated electron ef- 
fective mass in the oxides of metals with s 2 electronic con- 
figuration, fig. [H shows that the mass slightly decreases 
as the ionic radii of the cations from the same group or 
the symmetry of the same-cation oxide increases. How- 
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FIG. 3: Electronic band structure of single and multi-cation oxides, (a) Conduction band dispersion calculated in the afo-plane 
and (b) along the [0001] direction in InGaZnCU. (c) The conduction band unfolded nine times that corresponds to the number 
of [0001] layers in the unit cell, k=-|, c is the lattice parameter. The resulting single free-electron-like band is in accord with 
the stepless total density of states (d), in states/eV. (e) Conduction band of In203 is given for comparison, (f) Tight-binding 
conduction band (solid line) calculated for one-dimensional atomic chain depicted above the plot. Two types of metal atoms 
(red and purple spheres) alternate with oxygen atoms (blue spheres) and only the nearest-neighbor hopping (3 is assumed. To 
illustrate the effective mass averaging, cf. eq. (1), the conduction bands for the corresponding single-metal oxide chains (dashed 
lines) are aligned with (ei + £2)/2. 
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FIG. 4: The key electronic features of n-type TCO hosts. High energy dispersion of the conduction band (left) is essential for 
both high mobility of extra carriers due to their small effective mass and low optical absorption. For complete transparency in 
the visible range, the inter-band transitions E„ and E c should be >3.1 eV, while the intra-band transitions as well as plasma 
frequency should be <1.8 eV. Such a dispersed conduction band along with the above- visible optical transitions at the F point 
are found in a variety of oxides with s 2 -cation(s). The electron effective mass calculated for different oxide phases shows little 
dependence on the oxygen coordination and is isotropic {5=(m a + mb)/2m c — 1). The highlighted areas are to guide the eye; 
red circles represent currently known TCO hosts. 
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TABLE I: Net contributions from the states of the atoms that belong to the X-Ol or Y2-O2 layers (X=In or Sc, Y=G&, Zn, 
Al, Cd and/or Mg) to the conduction band wavefunction at V point, in per cent; and the electron effective masses m, in m e , 
calculated from the band structure of the layered oxides and the components of the electron effective-mass tensor, m a ,b and 
m z , calculated from eqs. (1) and (2) using the effective masses of the corresponding single-cation oxides. 



XY204 


Nx 


N01 


Ny 2 


N02 


m [100] 


ni[oio] 


m [001] 


mob 


m z 


InGaZn0 4 


23 


25 


29 


23 


0.23 


0.22 


0.20 


0.23 


0.23 


InAlCd0 4 


27 


27 


18 


28 


0.26 


0.25 


0.20 


0.27 0.27 


InGaMgCU 


27 


31 


21 


21 


0.27 


0.27 


0.24 


0.28 


0.29 


InAlMgCU 


33 


40 


12 


15 


0.32 


0.31 


0.35 


0.31 


0.34 


ScGaZnCU 


8 


19 


40 


33 


0.33 


0.33 


0.34 


0.33 


0.53 



ever, variations in the oxygen coordination as well as 
strong distortions in the metal-oxygen chains in different 
oxide phases lead to insignificant changes in the effec- 
tive mass. For example, for cubic (octahedral coordina- 
tion) and hexagonal (tetrahedral) ZnO or for corundum 
(distorted tetrahedral) and monoclinic /3-phase (both dis- 
torted tetrahedral and trigonal) Ga2C>3 the correspond- 
ing electron effective masses vary by 15% or 9%, respec- 
tively. The largest deviation in the effective mass values 
for various SiC>2 phases is 26%. Furthermore, the effec- 
tive mass remains isotropic for all phases of the s 2 -cation 
oxides - including silica ITQ-4 zeolite with large pore 
channels [24|]. These observations explain the success of 
amorphous transparent conducting oxides - in marked 
contrast to the amorphous Si where the directional inter- 
actions between the conduction p-orbitals lead to strong 
anisotropy of the transport properties which are sensitive 
to the orbital overlap and hence to the distortions in the 
atomic chains [14j . Finally, we note that the fact that the 
calculated as well as the observed [25| isotropic effective 
mass in rutile S11O2 where the direct overlap between Sn 
s-orbitals is possible (only) along the [001] direction, cor- 
roborate our conclusion that the s-s interactions do not 
govern the transport properties as discussed above. 

Effective mass averaging. Because the local symme- 
try (nearest neighbors), the cation-oxygen bond lengths 
and, hence, the s-p overlap are similar in the single and 
multi-cation oxides, the intrinsic transport properties in 
the layered materials should be related to those in the 
single-cation oxides. Moreover, due to the hybrid nature 
of the conduction states in the multicomponent oxides, 
the states of all cations should give the same order of 
magnitude contributions to the effective mass. Thus, we 
expect the later to be an "effective" average over the ef- 
fective masses of the corresponding single-cation oxides. 
Indeed, we find that in case of one-dimensional atomic 
chain where two types of metal atoms alternate with oxy- 
gen atoms, the effective mass averaging can be derived 
analytically, fig.[3£f). 

We formulate a simple approach which allows one to 
estimate the effective mass of the multicomponent ox- 



ides as follows. With proper doping, the Fermi energy 
is shifted up into the conduction band (Burstein-Moss 
shift). When the extra electrons propagate along the z 
direction, i.e., across the layers, the resulting resistivity 
is a sum of the resistivities of each layer. Therefore, the 
z component of the average effective-mass tensor can be 
found as: 

m z = (mi + m 2 + m 3 )/3, (1) 

where to 1,2, 3 are the effective masses of the corresponding 
single metal oxides - I112O3, Ga2C>3 and ZnO in the case 
of InGaZn0 4 . 

For the in-plane charge transport in the layered mate- 
rials, the effective- mass tensor components can be found 
in a parallel fashion. Note, that one needs to average the 
effective mass for the mixed GaZnC>2.5 layers: 



m a . b 3 \mi \{m 2 + m 3 ))' 

The resulting m a .b and m z are presented in table HI We 
find that the increase of the electron effective masses 
in the order InGaZn0 4 < InAlCd0 4 < InGaMg0 4 < 
InAlMg0 4 < ScGaZn0 4 is well reproduced by the above 
averaging. Moreover, the m Qj & and m z values nearly coin- 
cide with the corresponding effective masses of the multi- 
cation oxides with the exception of the Sc case. In fact, 
for the Sc-containing compounds the effective mass aver- 
aging is not legitimate due to the presence of the empty 
d-states of Sc near the conduction band edge. In SC2O3, 
the Sc d-states are located at 0.5 eV above the conduction 
band edge, and therefore give significant contributions to 
the effective mass which is about the mass of an electron. 
In the multi-cation oxide, the Sc d-states are found to be 
at ^2 eV above the conduction band edge, so that the 
resulting small effective mass, table HI is determined pri- 
marily by the s — p interactions - similar to the rest of 
the s 2 -cation oxides. 

The effective mass averaging procedure, eqs. (1) and 
(2), can be generalized for materials consisting of any 
number of layers, e.g., for InGa03(ZnO) m , TO=integer. 
Furthermore, since the intrinsic transport properties are 
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determined entirely by the local symmetry, i.e., the near- 
est neighbors, the effective mass averaging should apply 
to TCO's in amorphous state. In this case, one needs 
to average the components of the effective-mass tensor, 
m am orph={m a +m b +m z )/3. 

Importance of carrier generation. Upon doping of a 
TCO host material, the resulting conductivity depends 
not only on the effective mass but also on the carrier 
generation mechanism, carrier concentration and carrier 
relaxation time. Doping of a structurally anisotropic 
material may lead to non-uniform distribution of the 
carrier donors and, therefore, the isotropic behavior 
of the host may not be maintained as, for example, 
in oxygen-deficient /3-Ga2C>3 [13] • In multicomponent 
InGa03(ZnO) m , different valence states (In 3+ and Ga 3+ 
vs Zn 2+ ) and oxygen coordination (octahedral for In vs 
tetrahedral for Ga and Zn) are likely to result in pref- 
erential arrangement of aliovalent substitutional dopants 
or oxygen vacancies. Consequently, an anisotropic mo- 
bility should be expected in the layered materials due 
to the spatial separation of the carrier donors and the 
layers where the extra carriers are transfered efficiently, 
i.e., without charge scattering on the impurities. While 
targeted doping can help make either or both struc- 
turally distinct layers conducting, the amorphous com- 
plex oxides readily offer a way to maintain isotropic 
transport properties. Indeed, experimental observa- 
tions that the mobility and conductivity are indepen- 
dent of the large variations in the composition in amor- 
phous InGa03(ZnO)„ with n<4 [13] and that the effec- 
tive masses of amorphous and crystalline InGaZn04 are 
nearly the same [28| support our conclusions. 

For efficient doping of wide-bandgap oxides such as 
MgO, CaO, Si02 and AI2O3, novel carrier generation 
mechanisms should be sought. A non-traditional ap- 
proach has already yielded promising results in calcium 
aluminates [29I . I30I l3l| - a conceptually new class of 
transparent conductors Multicomponent oxides - 

such as those considered in this work, ordered ternary 
oxides [53 ] as well as their solid solutions and amorphous 
counterparts - represent an alternative way to utilize the 
abundant main-group elements such as Ca, Mg, Si and 
Al towards novel TCO hosts with a predictable effective 
mass and optical and transport properties controllable 
via the composition. 
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APPENDIX 

Theoretical methods First-principles full-potential linearized 
augmented plane wave method [33l . |33| with the local density ap- 
proximation is employed for electronic band structure investiga- 
tions of the XY2O4 compounds, X = In or Sc and Y = Ga, Zn, 
Al, Cd and/or Mg, and single-cation oxides. Cut-offs of the plane- 
wave basis (16.0 Ry) and potential representation (81.0 Ry), and 
expansion in terms of spherical harmonics with I < 8 inside the 
muffin-tin spheres were used. Summations over the Brillouin zone 
were carried out using 14 special k points in the irreducible wedge. 

Structure optimization. XY2O4 compounds have rhombohedral 
R3m crystal structure of YbFe2 O4 type [35l . l36fl . Indium (or scan- 
dium) atoms substitute Yb in 3(a) position, while both Y 3 + and 
Y 2 + atoms replace Fe in 6(c) position and are distributed ran- 
domly |37ll . Our total energy calculations for several structures in 
the (a,2a,c) supercell with various arrangements of the Y 2 + and 
Y 3 + atoms suggest that their arrangement is not ordered but ran- 
dom - in agreement with the experiment. We note here that the 
electronic band structure features are similar for different spatial 
distributions of the Y atoms for the reasons discussed in the paper. 

Since the valence state and ionic radii of Y 2 + and Y 3 + are dif- 
ferent, the site positions of these atoms as well as their oxygen sur- 
rounding should be different. Because the exact internal positions 
of atoms are unknown, we used those of the YbFe2 04 I'iTl as the 
starting values and then optimized each structure via the total en- 
ergy and atomic forces minimization. During the optimization, the 
lattice parameters were fixed at the experimental values |36l |38| . 
We find that the optimized cation-anion distances correlate with 
ionic radii of the cations. 

Single-cation oxides. For the single-cation oxides, the following 
phases have been calculated: Fm3ra for MgO, CaO, SrO, BaO 
and YbO; Ia3 for SC2O3 and Y2O3; Fm3m and P637TIC for ZnO; 
Fm3m for CdO; R3c for AI2O3; R3c and C2/m for Ga 2 3 ; Ia3, 
R3c and 72i3 for ln 2 3 ; P3 2 21, P2 1 /c, P6 3 /mmc, P4j2i2 and 



12 /m for SiC>2; P4i2i2 and P^/mnm for GeC>2; P^/mnm and 
Pbcn for Sn02- For each structure, the internal positions of all 
atoms have been optimized via the total energy and atomic forces 
minimization, while the lattice parameters were fixed at the exper- 
imental values. 

One-dimensional model of complex oxides. The effective mass 
averaging, cf. eq. (1), can be shown analytically using a one- 
dimensional model in the tight-binding approximation. To capture 
the key features of complex oxides, we consider a chain consisting 
of two types of metal atoms which alternate with oxygen atoms, 
fig. (3ff), and assume only the nearest-neighbor interactions given 
by the hopping integrals f3\ and /32 • The Hamiltonian of this model 
system is: 



H = V" \n,l)e t (n,l\ + V" \n',l')Pi(n,l\. 



(3) 



n,n' ,1,1' 



Here, I is the atom index in the unit cell, n enumerates the cells 
and n',1 1 in the second sum run over the nearest neighbors. For 
the bottom of the conduction band, the dispersion relation can be 
simplified to 



s(k). 



£1 + £2 



1 



If A , An 



[ka\ 



(4) 



if | £l -e a | < 2 



Here £0, £1 an d £2 are the atomic level 



energies of the oxygen and two types of metal atoms, respectively, 
and it is assumed that £0 < £1,2 and £1 ~ £2; A = \{e\ + £2) — £0 
and a is a half of the lattice parameter. Similar considerations 
for the chain consisting of only one type of metal atoms alternat- 
ing with oxygen atoms show that the quantity -|j represents the 

effective mass of the system. Therefore, eq. (4) represents the ef- 
fective mass averaging over those of the corresponding single-metal 
"oxide" chains, cf., fig. \3j[f) — in agreement with the results of our 
first-principles calculations. The following parameters were used to 
plot fig.EIf): £ = 1.00, £1 = 2.00, £2 = 2.05, /3i = 0.4 and ft> = 
0.5. 



